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The shape of the cross section of a dense fiber bundle is related to the symmetry of its molecular 
packing. However, this statement might be belied by type I collagen fibrils which have a rounded 
section of high symmetry while structural studies suggest that their molecules are assembled with a 
long range lateral order of lower symmetry. We examine how phyllotEixis, which is a non conventional 
crystallographic solution to packing efficiency in situations of high radial symmetry, can establish 
a link between those two apparently conflicting points. The lateral order imposed by the algorithm 
of phyllotaxis, which implies an enlargement of the notion of long range lateral order beyond that 
used for classical crystals, provides a basis for a new analysis of the experimental data. 
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FIG. 1. X-ray scatterings along the direction normal to the fibril axis, equatorial traces, obtained from native tendons by 
Doucet et al.(2011) (a), Hulmes et al.(1981) (b) and MacBride et al.(1997) (c). We present here only raw equatorial traces of the 
scattered intensities. They have not been submitted to any computational treatment aiming at suppressing their background 
scatter in order to isolate the small discrete peaks merging out from it in b) and c). 

INTRODUCTION 

The organization of long molecules of biological origin in dense fiber bundles can present various degrees of order 
according to the nature of the material and its environment. For instance DNA molecules can be hexagonally packed 
in large toroidal condensates with hexagonal cross section [J Q while type I collagen molecules can show a high level 
of disorder in fibrils with circular cross section Q. Those two extreme examples illustrate the relation which can 
be expected between the local symmetry, as imposed by the packing of the molecules, and that of the aggregate, as 
shown by the cross section of the bundle. As a matter of fact, if the molecules are packed with a long range lateral 
crystalline order the cross section of their bundle is limited by the polygonal trace of reticular planes with lowest 
interfacial energy, if not this energy is minimized by a circular cross section. However, this clearcut point of view has 
to be reconsidered as type I collagen fibrils with circular cross section are not always as disordered as those mentioned 
above. X-ray experiments indeed suggest that such fibrils might present some long range lateral order as shown by 
the equatorial traces of the raw scattered intensity obtained from such fibrils where they are compared with that 
of a disordered fibril Q in figure [1] 

These traces all present an important diffuse scattering which, after having decreased out of the beamstop, attains 
a broad maximum for fc w 5 nm~^. On trace (a) this diffuse scattering is smooth while on traces (b, c) small slightly 
broadened peaks merge out from it in the 1 < A: < 5 nm"-'^ range as shown by vertical arrows. The maximum at 
/c « 5 nm~^ and the global shape of the diffuse scattering correspond well to a highly disordered organization with 
a mean intermolecular distance of 1.3 nm, but the peaks observed for 1 < k < 5 nm~^ suggest the presence of some 
kind of order with repeat distances larger than the intermolecular one. Thus these observations reveal a complex 
material in which disorder, the diffuse scattering, and some kind of order, the discrete bumps, can coexist without the 
reticular planes of the latter manifesting themselves in a facetted cross section with low symmetry. This complexity 
was indeed pointed out for the first time forty five years ago0| and is regularly addressed since then 

CONCILIATING ORDER AND DISORDER IN TYPE I COLLAGEN FIBRILS 

The data are usually analyzed considering that the formation of type I collagen fibrils proceeds along two steps 
as shown in figure [21 The segments of triple helices are densely organized in an overlap region but in a gap region 
one segment of triple helix out of five is replaced by a vacancy and the distribution of the vacancies form periodic 
layers of Hodge-Petruska staggering. Reviews and analysis of the large number of models proposed to interpret the 
experimental data on the basis of this description can be found in0, [l^. Those models are of very different natures, 
going from a purely crystalline organization to a liquid-like one, including multiple-start spirals and concentric ring 
models. A rather satisfying agreement is obtained considering that : 

- the vacancies and the triple helix sections at the gap/overlap interface are organized according to the x,y plane of a 
triclinic crystal to give account of the discrete scatterings, 

- the gap and overlap regions are disordered, eventually more in the first, to give account of the background scatter, 

- the whole fibril is described as a polycrystalline assembly of monocrystalline grains with triangular basis separated 
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1.3 nm 




FIG. 2. Left-handed collagen molecules (a) assemble in right-handed triple helices (b,c) represented by rods (d) which in turn 
are proposed to be assembled with a regular shift (e), the so-called Hodge-Petruska staggering, in order to create regions of 
gap "g " and overlap " o " to give account of the longitudinal striations of the fibrils. 



by radial grain boundaries in order to approach a circular section by multiplying the number of monocrystalline facets 
at the periphery of the poly crystal. 

We recently proposed to look at this structural problem along a quite different way inspired by the iterative process 
of phyllotaxis, a non conventional crystallographic solution to packing efficiency in situations of high radial symmetry 
[llLll2l|. This process ensures to each of the molecules the most homogeneous and isotropic local environment possible, 
at the price of local metric and topological disorders which interact to build a long range self-similar organization 
with high radial symmetry. In 11 1 we considered parallel triple helices in fibrils of small section, represented by a 
phyllotactic pattern of 2500 points, and it appeared that the calculated X-ray scatterings presented features which 
could contribute to that observed along the direction normal to the fibril axis in the 1 < fc < 5 nm^^ range. 

We develop here this preliminary approach considering phyllotactic patterns with larger numbers of points, up to 
15704 approaching the most common sizes of the fibrils, and making full use of our recent study of the symmetries 
presented by such patterns [l2|. This development shows that a long range order can be compatible with a circular 
cross section, an order of course different from that of classical crystallography, and that the scattered intensity 
expected from it can be compared with the localized scatterings observed in figure [T] We also introduce random 
displacements around the positions of these points to simulate some complementary disorder and give a better account 
of the background diffuse scattering observed in the same figure. Moreover, this phyllotactic order imposes an intrinsic 
structural complexity which suggests a new reading of the experimental data, particularly their eventual variability 
with the size of the cross section. 



The fan shape of the scattered intensities out of the equatorial direction suggests the presence of a twist in the assembly of the triple 
helices, it is however weak and we ignore at this stage of study. 
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PHYLLOTAXIS 



The lateral organization of parallel fibers can be represented by their projection onto a plane normal to the axis 
of the bundle. We build the phyllotactic organization in this plane. In the case of type I collagen two situations are 
to be considered corresponding respectively to the overlap regions, in which all points are occupied, and to the gap 
regions, in which one point over five is a vacancy. 



Phyllotactic pattern for an overlap region 

A phyllotactic organization of points indexed by integer s is built by an algorithm such that the position of point 
s is given by its polar coordinates r = a^/s and 9 = 2ttXs that is r = a^/O/ -\/27rA which is the equation of a Fermat 
spiral here called the generative spiral. The area of the circle of radius r which contains s points is tto^s so that 
the area per point has the value -ko? , indeed it oscillates close to this value for small s then converges towards it 12 1. 



The most homogeneous and isotropic environment, or the best packing efficiency in radial symmetry, is obtained with 
A = 1/r where t is the irrational golden ratio (1 + ■\/5)/2 |13l4l5| . A sector of a phyllotactic pattern for N = 15704 
points with their Voronoi cells is shown in figure |3l 

In such a pattern, pentagons and heptagons are topological defects distributed among the hexagons. They appear 
concentrated in narrow circular rings with constant width which separate large rings of hexagons whose width increases 
as one moves from the core towards the periphery. In the narrow rings, pentagons and heptagons are associated in 
dipoles separated by hexagons whose shape is close to that of a square with two corners cut. The rings of dipoles 
are indeed grain boundaries separating hexagonal grain^. The underlying arithmetic of this organization is the 
Fibonnacci series, = fu-i + fu-2 from /o = and /i = 1, which makes the organization self-similar as it is 
invariant by a change of scale r" as summarized in table msee[l3). This is also made apparent considering the 
evolutions of the distances between first neighbors shown in figure 21 



Phyllotactic pattern for a gap region 

The phyllotactic pattern of a gap is similar to that of an overlap except for the fact that one point out of five must 
be replaced by a vacancy. The very nature of the iterative process leads to do this replacement all along the generative 
spiral. Then, the vacancies appear distributed on new spirals in a manner respecting the periodicity required by the 
staggering of the triple helices, as shown in figure[S] On these new spirals, the vacancies can be first neighbors or not. 
They are first neighbors in the hexagonal grains having one 5s multiple of five but not in all others. This organization 
can be built from table |T] starting from one point s in a hexagonal grain and using the 6s of this grain to number its 
surrounding sites step by step. Two examples are drawn in figure |6] based on regular hexagons for a point s in grain 
u — 10, where \5s\ — 89, 55, 34 with one 5s multiple of 5, and a point s in grain 12 where \5s\ = 233, 144, 89 without 
5s multiple of 5. If the starting point s of figure [6] is considered as a vacancy and its number s is a multiple of 5, the 
vacancies are first neighbors in grain 10, second neighbors in grain 12 and build lattices whose elementary cells are 
shown in figure [7] in these drawing based on regular hexagons. 

In each case, (a) or (b) , the staggering of the triple helices shown in figure [2^ can be lined up along straight lines 
joining first neighbors in agreement with what is observed along the spirals of the parastichies in the pattern of 
figure m Of course, in a grain of the phyllotactic pattern, the Voronoi cells are hexagonal in the topological sense 
only and the elementary cells (a) and (b) drawn in figure [7] suffer deformations when moving in between the grain 
boundaries limiting the grain. The notion of elementary cell is no longer that used in classical crystallography. It 
has lost its metric content, distances and angles defining the positions of the sites, to keep only its topological part, 
neighboring relations. 

As shown in table |lll those elementary cells (a) or (b) are alternatively localized in concentric domains of increasing 
sizes corresponding to groups of three consecutive grains, domains A, and groups of two consecutive grains, domains 
B. As in domains A one number 5s of parastichies is a multiple of 5 the simple elementary cell (a) can be periodically 
packed in the circular grains and the vacancies are first neighbors. This is no longer possible in domains B where the 
numbers of parastichies 5s are not multiple of 5 and impose the more complex elementary cell (b) . 



^ The grain boundaries are here circular while they are radial in the polycrystalline model of Hulmes et al. (1995) a dense organization 
with radial symmetry can not be built without defects. 
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FIG. 3. A quadrant of a set of 15704 points organized according to the algorithm of phyllotaxis with the golden ratio. The radius 
of this quadrant is about 80nm. Each point is surrounded by its polygonal Voronoi cell whose number of sides corresponds 
to the number of first neighbors around this point. Blue cells are pentagons, orange cells hexagons and green cells heptagons, 
those polygons are not necessarily regular. The three spirals joining the first neighbors of the points are called parastichies. 
In the case of a point in a hexagonal Voronoi cell three spirals, called parastichies join the first neighbors of the point. These 
parastichies are Fermat spirals, which are not to be confused with the Archimedean spirals of the spiral models described in 
Hulmes et al. (1995). 



In a domain A such as that of grains 9, 10, 11 where Ss = 55 there are 55/5 — 11 spirals of first neighbor vacancies. 
These spirals remain in the domain B of grains 12, 13, but vacancies are now second neighbors and the elementary 
cell is of type (b). In this domain the 11 spirals are aligned along edges of the elementary cells. In the next domain 
A of grains 14, 15, 16 the 610/5 = 122 spirals would appear formed with first neighbor vacancies. Nevertheless they 
are discernable before in grains 12, 13 formed with second neighbors. These spirals are aligned along edges of the (b) 
cells. So in these grains one family of cell edges form 11 spirals and the other family form 122 spirals. 



PHYLLOTAXIS VS CLASSICAL CRYSTALLOGRAPHY 



In a phyllotactic pattern, most of the points have six first neighbors, as expected for a dense system, but those 
first neighbors can not draw regular hexagons. They are distorted like in a flat 2D crystal bent in its plane so that 
fluctuations of distances and angles introduce a metric disorder. Moreover, a few points have flve or seven flrst 
neighbors which are defects, disclinations, introducing a topological disorder. Those two disorders interact so as to 
build a concentric organization of grains of points with six first neighbors separated by grain boundaries in which the 
topological defects arc associated in dipolcs equivalent to dislocations. Those dislocations arc distributed according 
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TABLE I. Cell types in the successive rings. The number s is that of each point on the generative spiral, the first neighbors of 
point s have numbers s + Ss where Ss are Fibonnacci numbers, all the same for a given ring except in the core of the pattern. 
The large hexagonal rings or grains are labeled by the rank u of the Fibonnacci number /„ corresponding to the medium value 
of 6s in a ring, the Ss are also the numbers of parastichies of each type in the ring. The narrow rings or grain boundaries are 
marked by ||. 



u 


Cell type 


number of cells 


s from 


to 


neighbor separations 5s 




pentagon 


2 





1 


1,2,3,4,5 or -1,2,3,5,8 




hexagon 


1 


2 


2 


-2,2,3,5,8,13 




heptagon 


3 


3 


5 


(-3,-2,2) or (-4,-3,-2) or (-5,-3,-2) ,3,5,8,13 




hexagon 


1 


6 


6 


-5,-3,3,5,8,13 




pentagon 


2 


7 


8 


-5,-3,5,8,13 




hexagon 


1 


9 


9 


-8,-5,-3,5,8,13 


II 


hexagon 


5 


10 


14 


-8,-5,5,8,13,21 


II 


heptagon 


3 


15 


17 


-13,-8,-5,5,8,13,21 


II 


hexagon 


5 


18 


22 


-13,-8,-5,8,13,21 


II 


pentagon 


8 


23 


30 


-13,-8,8,13,21 


7 


hexagon 


2 


31 


32 


-21,-13,-8,8,13,21 


II 


heptagon 


13 


33 


45 


-21,-13,-8,8,13,21,34 


II 


hexagon 


8 


46 


53 


-21,-13,-8,13,21,34 


II 


pentagon 


13 


54 


66 


-21,-13,13,21,34 


8 


hexagon 


34 


67 


100 


-34,-21,-13,13,21,34 


II 


heptagon 


21 


101 


121 


-34,-21,-13,13,21,34,55 


II 


hexagon 


13 


122 


134 


-34,-21,-13,21,34,55 


II 


pentagon 


21 


135 


155 


-34,-21,21,34,55 


9 


hexagon 


134 


156 


289 


-55,-34,-21,21,34,55 


II 


heptagon 


34 


290 


323 


-55,-34,-21,21,34,55,89 


II 


hexagon 


21 


324 


344 


-55,-34,-21,34,55,89 


II 


pentagon 


34 


345 


378 


-55,-34,34,55,89 


10 


hexagon 


422 


379 


800 


-89,-55,-34,34,55,89 


II 


heptagon 


55 


801 


855 


-89,-55,-34,34,55,89,144 


II 


hexagon 


34 


856 


889 


-89,-55,-34,55,89,144 


II 


pentagon 


55 


890 


944 


-89,-55,55,89,144 


11 


hexagon 


1221 


945 


2165 


-144,-89,-55,55,89,144 


II 


heptagon 


89 


2166 


2254 


-144,-89,-55,55,89,144,233 


II 


hexagon 


55 


2255 


2309 


-144,-89,-55,89,144,233 


II 


pentagon 


89 


2310 


2398 


-144,-89,89,144,233 


12 


hexagon 


3384 


2399 


5782 


-233,-144,-89,89,144,233 


II 


heptagon 


144 


5783 


5926 


-233,-144,-89,89,144,233,377 


II 


hexagon 


89 


5927 


6015 


-233,-144,-89,144,233,377 


II 


pentagon 


144 


6016 


6159 


-233,-144,144,233,377 


13 


hexagon 


9167 


6160 


15326 


-377,-233,-144,144,233,377 



to an inflation/deflation quasicristalline rule along the grain boundaries whose radii vary as the Fibonacci sequence. 
This organization is therefore not invariant by the translation and rotation operations of classical crystallography, 
but is self-similar. The grain boundaries virhich, in ordinary materials, are extrinsic defects imposed by deformations 
of the material arc here intrinsic. In the grains, the six first neighbor points of each point are aligned along the three 
families parastichies which can be seen as equivalent to the reticular planes of classical crystallography. Crossing a 
grain boundary, from one grain to the next, the dislocations present in it introduce the new parastichies, hence the 
new points, needed to maintain the density constant (see 12| for remarks in this paragraph). 

The staggering of the triple helices introduces vacancies which are organized along new spirals in each gap where 
they can be either first neighbors or not. Each of these two distributions can be said ordered, but once again not 
in the sense of classical crystallography, as this order concerns only the neighboring relations of one vacancy among 
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FIG. 4. Distances between first neighbor points with the parameter a — 0.71 leading to a mean distance close to 1.3 nm as in 
collagen. The three colors correspond to the three parastichies, the upper and lower crossings on the same verticals correspond 
to grain boundaries and the intermediate ones to the cores of hexagonal grains. All distances are confined in the same domain 
whatever the size of the structure from 1.18 nm to 1.77 nm. 



TABLE 11. Grains from u = 8 to ti = 21, 5s in those grains (multiple of 5 are in bold face), domains A and B occupied by 
elementary cells (a) and (b), number of points in the domains. 



Hexagonal 


s from 


to 


Ss 


Concentric 


inner 


external 


Number of 


grains 








domains 


radius (nm) 


radius 


points in domain 


9 


156 


289 


55,34,21 


Al 


8.5 






10 


379 


800 


89,55,34 


A 1 






2010 


11 


945 


2165 


144,89,55 


AJ 




31.7 




12 


2399 


5782 


233,144,89 


Bl 


33.43 






13 


6160 


15326 


377,233,144 


BJ 




84.25 


12928 


14 


15937 


40428 


610,377,33 


Al 


86.16 






15 


41416 


106335 


987,610,377 


A 1 






263250 


16 


107933 


279186 


1597,987,610 


AJ 




360.6 




17 


281771 


732208 


2584,1597,987 


Bl 


362.3 






18 


736390 


1919036 


4181,2584,1587 


BJ 




945.5 


1637266 


19 


1925822 


5027483 


6765,4184,2584 


AJ 








20 


5038430 


13167594 


10946,6765,4184 


A 1 








21 


13185306 




17711,10946,6765 


AJ 









four cross sections of triple helices. In one case, elementary cell (b), the disposition of the four cross sections of triple 
helices around one vacancy is very close to that existing in the x,y plane of the triclinic cell proposed in some previous 
models However the fluctuations of distances and angles characteristic of phyllotaxis can not let expect Bragg 

peaks in the X-ray diagrams. Finally, those two distributions of vacancies are themselves organized in alternated 
concentric domains whose sizes increase with the radius of the pattern, in a manner related to the Fibonacci series. 
This is an intrinsic structural heterogeneity which has no equivalent in classical crystallography. 

SCATTERED INTENSITIES EXPECTED FROM PHYLLOTACTIC PATTERNS 

The main information concerning the lateral organization of collagen fibrils being provided by X-rays scattering 
experiments, we calculate the scatterings expected from the above phyllotactic templates for comparison. As the 




FIG. 6. Numberings of the points surrounding a point s in grain 10 with one 5s multiple of 5 (a) and grain 12 without 5s 
multiple of 5 (b), for simpleness the Voronoi cells of the grains are here assumed to be regular hexagons and the parastiches 
are therefore straight lines. 
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a) b) 



FIG. 7. Elementary cells determined by the vacancies, open circles, in the triangular network of points of figure [S] for grain 10 
with one Ss multiple of 5 (a) and grain 12 without Ss multiple of 5 (b). The multiplicity of both cells in the hexagonal lattice 
is 5. 



FIG. 8. Intensities scattered by overlaps, taking account of a form factor corresponding to a molecular diameter of 1 nm, for 
N=2500 (a), N=6500 (b) and N=15704 (c). 




triple helices are normal to the patterns, the equatorial traces of the intensities scattered by the cylindrical overlaps 
and gaps containing them are those of their respective patterns. The scattered intensity expected from a distribution 
of points is the square of the amplitude of its Fourier transform, but, as the points are representative of triple helices 
having a diameter of 1 nm, the points are to be changed into disks and the intensity decreases as the form factor of 
the disk. The results are shown in figures |8] and |9] for the overlap and gap patterns respectively, in both cases for 
patterns of increasing sizes. 

The intensity scattered by overlaps show that, whatever the size of the pattern, the scattering is concentrated in the 
two regions fc < 1 nm^^ and A: w 5 nm~^. The scattering in the small wave vectors region corresponds to long range 
fluctuations such as the size of the pattern and the circles of dipoles. The scattering around 5 nm^^ is to be related 
to the local organization of the triple helices. In the intensity scattered by the gaps, the features of the scatterings 
produced by the overlaps are still present but new ones appear at A: ~ 1 nm^^ and in the range 1 < k < 5 nm^^, 
the latter, hardly discernable for N=2500, become more important for N ~ 6500 then quite obvious for N = 15704. 
These are not strictly Bragg peaks, but peaks with a finite width and broadened feet which contribute to a raising of 
the baseline. 



FIG. 9. Intensities scattered by gaps, taking account of a form factor corresponding to a molecular diameter of 1 nm, for 
N=2500 (a), N=6500 (b) and N= 15704 (c). 




» 1 2 3 4 5 «1,„™-1 1 2 3 4 5 6 




1 2 3 4 5 6 
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TABLE III. Positions measured on raw equatorial traces (A) Hulmes 1981,1995, Mc Bride 1997 , on traces from which the 
background scatter has been subtracted (B) Wess 1998, Orgel 2001 and on the calculated trace (C) of figure 9c. 





peaks (nm ^) 


A 


ref. [4] 




1.65 


2.5 


3.3 


3.59 


4.6 


4.98 






ref. [5] 




1.62 


2.65 


3.27 


3.54 


4.5 


4.79 






ref. [6] 




1.64 


2.5 


3.28 


3.57 


4.58 


4.92 




B 


ref. [101 




1.4 


2.1;2.3 


3.19 


3.47 


4.67 


5.17 


(5.7) 




ref. [161 




1.3 


2.2;2.3 


3.2 


3.47 


4.67 


5.04 


(5.7) 


C 


Phyllotaxis 


1.05 


1.6 


2.4 


3.2 


3.4 


4.5 


5 


(5.6) 



FIG. 10. Scattered intensity obtained from a pattern having the size of a pattern with N=15704 points but in which the only 
scattering centers are the points corresponding to the vacancies. 




The scatterings observed at /c ~ 1 nm^^, one fifth of the position of that observed at fc ~ 5 nm~^, is that expected 
from elementary cells (a) and those observed in the 1 < fc < 5 nm^^ range are those expected from elementary cells 
(b). The evolutions of these scatterings with the sizes of the patterns can be understood from table ITTl The patterns 
with TV = 2500 is dominated by grains 9 to 11, i.e. a domain A with the elementary cell (a), increasing the size of 
the pattern to iV = 6500 then N ~ 15704 adds grains 12 then 13, both included in a domain B with elementary cell 
(b) whose contribution increases accordingly. If N was increased beyond 15704 a domain A would be added whose 
contribution would become rapidly dominant as it contains about 20 times more points than the two first domains. 

COMPARISON WITH EXPERIMENTS 
X-rays data 

Raw equatorial traces obtained with native samples have been already shown in figure [T] and described in the 
introduction. The positions of the peaks merging out from the the background scatter in figures lb,c reasonably agree 
with those obtained for a gap represented by a phyllotactic pattern of N=15704 points, as shown in table Hill Such 
an agreement is not surprising as the distribution of the vacancies in the domains B of the gaps is, at least as far 
as its topology is concerned, close to that in the crystalline cells proposed in previous works. The widths of these 
discrete scatterings results from the metric distortions of the Voronoi cells and their distribution within the frame of 
the self-similar reproduction. The scattering of domain A at fc = 1 nm^^ has not been mentioned up to now. A trace 
of it might be found in figure [IJi from ref. [Sl but may be lost in the intense small angle scattering observed around 
the beamstops in the other references. 

When heavy atoms are incorporated in the sample in order to enhance its contrast the general shape of the scattered 
intensity changes with respect to that of native samples Essentially the very intense scattering at fc « 5 nm~^ 
decreases and that at fc sa 3 nm~^ increases. This change might result from a change of structure induced by the 
incorporation process, but it also might be that the heavy atoms are predominantly localized in the vacancies of the 
gaps enhancing the weight of the scattering expected from those vacancies as shown on figure [TOl 

But, in the experimental data, those discrete scatterings contribute for no more than 5 to 10 percent of the total 
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FIG. 11. Scattered Intensities by N=15704 sites without disorder (a), with an increasing disorder (b,c). The histogram of 
radial displacements is indicated (azimuthal displacements are similar). These calculated intensities take account of an equal 
contribution from overlaps and gaps. 



2 4 





* nm 




intensity [8[ which is dominated by a background scatter bearing evidence for an important level of disorder in the 
sample. Certainly, a phyllotactic organization contains an intrinsic metric disorder associated to the deformations of 
the Voronoi cells, but as shown by the calculated scattered intensities presented in figures [5] and [HI these fluctuations 
are not sufficient to give account of the intense background scatter observed. Indeed, phyllotaxis is but a geometrical 
template and a disorder must be introduced. Such a disorder can be simulated d making the position of each point 
s fluctuates randomly 17, 3] with an amplitude characterized by the histograms of figure [TT] where the evolution of 
the calculated spectra as this disorder increases is shown. 

A randomly distributed disorder does not change the shape of the scattering, it just change the amplitude and adds 
a constant background. It is analogous to what happens when such disorder is introduced in a crystal: the width 
of Bragg peaks remain very small and a Laue continuous background appears, but if correlations are introduced in 
the disorder they lead to a broadening of localized peaks. To introduce short range correlations we have convoluted 
displacement with a Gaussian function having a width close to the distance between second neighbors. This method 
introduces short range correlations in the disorder between first and second neighbor points but not long distance 
ones, nevertheless the broadlines of the experimental spectra are rather well reproduced: small discrete bumps are 
superposed on an background increasing with k and culminating at fc w 5 nm~^. The positions of the small discrete 
peaks are not affected by the disorder, only their width, this holds to the fact that the displacements on short distances 
are small compared to the characteristic distances between rows of vacancies. 



Electron microscopy 



Micrographs of ultrathin sections of rat tail tendon prepared to provide optimum preservation of the native fibril 
structure with sufficient contrast for electron microscopy give an impression of tangentially oriented lines of density 
which may be interpreted as concentric circles or spirals The study by optical diffraction with an equivalent 
aperture of 120 — 180 nm of micrographs of fibril sections with large radii of 210 nm measures a separation period of 
4.6 nm and indicates that this periodicity is detected predominantly towards the edge rather than the center of the 
fibrils. This can be put in correspondence with the heterogeneous organization of a phyllotactic pattern of similar size 
as deduced from table |TT1 In such a pattern a peripheral ring of thickness 120 nm made of a domain A surrounds a 
central core of radius 80 nm made of a domain B predominantly. In the peripheral domain A first neighbor vacancies 
outline 122 spirals distant by 3 to 5 nm and arriving at a small angle on the limit of the pattern while in the central 
domain B second neighbor vacancies outline a more complex combination of spirals. 



CONCLUSION 



Electron microscopy studies show that type I collagen fibrils have a rounded cross section while some X-rays 
scattering studies suggest that their components, the so-called triple helices, can be assembled with a long range 
lateral order. This order cannot be crystalline as the topological constraint of circular symmetry would not be 
compatible with its long range propagation. The most homogeneous and isotropic dense organization compatible 



^ We do not establish a distinction between gaps and overlaps at this stage although the disorder might be more important in the first 
because of the vacancies. 
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with circular symmetry is indeed the self similar one issued from the iterative process of phyllotaxis. We considered 
here that the circular symmetry is imposed to the assembly of triple helices, for instance by the interfacial tension 
of the collagen fibril, and that this assembly adapts itself to this constraint in order to achieve the best packing 
efficiency typical of phyllotaxis. This adaptation while the fibril grows requires a certain internal mobility which is 
indeed suggested by NMR studies (lit- 

When a long range lateral order is detected by X-rays scattering studies, it is characterized by distances larger 
than the mean distance between triple helices. This is most likely to be associated with the fact that, according 
to the Hodge-Petruska staggering, one section of triple helix out of five in the gaps of the fibrils is replaced by a 
vacancy. Doing this replacement all along the generative spiral of phyllotaxis leads to a distribution of vacancies in 
the phyllotactic pattern of the gaps such that the staggering periodicity is established along parastichies of the pattern 
following an organization in concentric alternated domains presenting two different structures. 

The scattered intensity expected from one of these structures is made of discrete scatterings whose positions compare 
well with those observed. The distribution of vacancies is indeed, at least as far as its topology is concerned, close 
to that in the crystalline cells proposed in previous works to analyze the data. However, here, although narrow in 
their most intense part, these scatterings are not Bragg peaks, they have a more complex shape due to the metric 
fluctuations inherent to phyllotaxis. The presence of domains with another structure has not been mentioned up to 
now. except may be suggested by electron microscopy observations. The characteristic component of the intensity 
scattered by this structure should be observed in the small angle region and might contribute to the very intense small 
angle scattering observed close to the beamstop in most of the observations. But also, as the widths of these domains 
increase rapidly when moving away from the center, the scattering of a fibril should be dominated by that of one or 
the other of the two domains according to its size. This could underlie the differences between the spectra shown in 
figure [2 Unfortunately, scattering studies do not give a direct access to the size of the fibrils and the samples are 
tendons which often contain fibrils of different sizes. 

Finally, phyllotaxis is but a zero temperature model, some thermal disorder must be introduced in the gaps as well 
as in the overlaps. Our simulation, although limited at local random displacements without long range correlations, 
leads to a rather good reproduction of the broadlines of the observed spectra. As those displacements are small 
compared to the characteristic distances between rows of vacancies, only the width of the scatterings expected from 
them is affected, not their position. This broadening increases with the disorder making these discrete scatterings less 
discernable, an evolution which might also contribute to the differences between the spectra shown in figure [TJ 
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